Using photoemission spectroscopy to probe a strongly interacting Fermi gas 
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Ultracold atom gases provide model systems 
in which many-body quantum physics phenom- 
ena can be studied. Recent experiments on 
Fermi gases have realized a phase transition to 
a Fermi superfluid state with strong interparti- 
cle interactions [3, 0, 0, H, H, @, 0]. This sys- 
tem is a realization of the BCS-BEC crossover 
connecting the physics of BCS superconductiv- 
ity and that of Bose-Einstein condensation (BEC) 
While many aspects of this system 
have been investigated, it has not yet been possi- 
ble to measure the single-particle excitation spec- 
trum, which is a fundamental property directly 
predicted by many-body theories. Here we show 
that the single-particle spectral function of the 
strongly interacting Fermi gas at T ss T c is dra- 
matically altered in a way that is consistent with 
a large pairing gap. We use photoemission spec- 
troscopy to directly probe the elementary exci- 
tations and energy dispersion in the Fermi gas 
of atoms. In these photoemission experiments, 
an rf photon ejects an atom from our strongly 
interacting system via a spin-flip transition to a 
weakly interacting state. We measure the occu- 
pied single-particle density of states for an ul- 
tracold Fermi gas of 40 K atoms at the cusp of 
the BCS-BEC crossover and on the BEC side of 
the crossover, and compare these results to that 
for a nearly ideal Fermi gas. Our results probe 
the many-body physics in a way that could be 
compared to data for high-Tc superconductors 
This new measurement technique for ultra- 
cold atom gases, like photoemission spectroscopy 
for electronic materials, directly probes low en- 
ergy excitations and thus can reveal excitation 
gaps and/or pseudogaps. Furthermore, this tech- 
nique can provide an analog to angle-resolved 
photoemission spectroscopy (ARPES) for prob- 
ing anisotropic systems, such as atoms in optical 
lattice potentials. 

As interacting quantum systems with highly tunable 
parameters and well understood two-body interactions, 
ultracold atom gases provide model systems in which to 
test condensed matter theories. A challenge for experi- 
menters is to find ways to probe these atom gases that 
relate directly to condensed matter ideas and enable sen- 
sitive searches for new phenomena that can advance our 
understanding of strongly correlated systems. At a very 
basic level, the effect of interactions is a modification of 




FIG. 1: Photoemission spectroscopy for ultracold 
atom gases, a In electron PES, one measures the energy of 
electrons emitted from solids, liquids, or gases by the photo- 
electric effect. Using energy conservation, the original energy 
of the electrons in the substance can be determined. Simi- 
larly, in photoemission spectroscopy for atoms, an rf photon 
with energy, hv, transfers atoms into a weakly interacting spin 
state, b The rf photon drives a vertical transition where the 
momentum hk is essentially unchanged. By measuring the en- 
ergy and momentum of the out-coupled atoms (upper curve) 
we can determine the quasiparticle excitations and their dis- 
persion relation (lower curve). Here cf> is the Zeeman energy 
difference between the two different spin states of the atom. 



the single-particle states. As interactions are increased, 
the single-particle eigenstates of the non-interacting case 
become quasi-particles and phase transitions manifest 
themselves as qualitative changes to the excitation spec- 
trum fl^ |. such as the appearance of energy gaps. The 
single-particle excitation spectrum can be predicted by 
many-body theory and is a fundamental property of any 
interacting system. 

For electronic systems, photoemission spectroscopy 
(PES) provides a powerful technique to probe the occu- 



pied single-particle states [13|. In a typical PES experi 



ment, electrons are ejected from a substance through the 
photoelectric effect, see Fig. QJi. The photoelectrons are 
collected, energy and momentum resolved, and counted 
to give a spectrum of intensity as a function of the mea- 
sured kinetic energy, = h 2 k 2 /2m. Here, % = h/2ir, 
where h is Planck's constant, and m is the particle mass. 
By conservation of energy, one can determine the energy 
of the original single-particle state, E s , using 



E s = £fe + <p ~ hv. 



(1) 



Here, hv is the photon energy, <f> is the work function 
of the surface, and Ep — E s is often referred to as the 
binding energy (l3j . 
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For ultracold atom gases, radio-frequency (rf) spec- 
troscopy has been used to probe a strongly interacting 
Fermi gas 0, 0, EI EI EE EE 0- In a typical ex- 
periment, a pulse of rf drives atoms into an unoccupied 
Zeeman spin state, where they are counted to yield a 
spectrum of counts versus rf frequency. To date, the rf 
out-coupled atoms have not been energy or momentum 
resolved. However, analogous to electron PES, the mo- 
mentum of the rf photon is negligible compared to the 
typical momentum of the atoms and therefore the mo- 
menta of the out-coupled atoms are characteristic of the 
original atom states. Eqn. [1] applies to photoemission 
spectroscopy of atom gases, by means of momentum- 
resolved rf spectroscopy, if one simply replaces the work 
function cj) with the Zeeman energy splitting, see Fig. Q]d. 
The extension of photoemission spectroscopy from con- 
densed matter to cold Fermi gases was discussed by Dao 
et al. EH. 

In this paper, we use photoemission spectroscopy, by 
means of momentum-resolved rf spectroscopy, to probe 
an ultracold gas of fcrmionic 40 K atoms. Similar to PES 
in solids, this measurement probes the single-particle 
spectral function, which is directly related to the single- 
particle Green's function predicted by many-body theo- 
ries 13(. We use this new technique to probe the Fermi 
gas near a magnetic-field Fano-Feshbach resonance where 
one can tune strong atom-atom interactions to realize a 
Fermi superfluid in the region of the BCS-BEC crossover 

Our Fermi gas consists of 3 x 10 5 40 K atoms in a mix- 
ture of two spin-states. The gas is confined in an optical 
dipole trap and evaporatively cooled to T/Tp = 0.18, 
where T is the temperature, Tp is the Fermi tempera- 
ture as defined by Tp = Ep/ks, and fcs is Boltzmann's 
constant. The Fermi energy, Ep — h ■ (9.4 ± 0.5 kHz), 
is determined from a measurement of the peak density 
of the trapped gas. For the photoemission spectroscopy, 
we apply an rf pulse that couples atoms in one of the 
two spin states to an unoccupied third spin state. There 
are two essential requirements for determining the ex- 
citation spectrum. The first is that the interaction en- 
ergy is sufficiently small that = h 2 k 2 /2m holds and 
the data are not subject to complicated final-state ef- 
fects [H ED, 0, El [ij, El El ■ The second requirement 
is that collisions do not scramble the energy and mo- 
mentum information carried by the out-coupled atoms. 
In previous rf spectroscopy measurements both of these 

requirements were not satisfied S El EI EI El- in 



our K gas, however, the interaction energy for the out- 
coupled atoms is approximately 640 Hz, which is much 
smaller than Ep. Furthermore, the mean- free path for 
the out-coupled atoms is much larger than the size of the 
gas: -1- « 6 Rf, where a is collision cross section, n is 
the average density, and Rp is the Fermi radius of the 
non-interacting gas. 

To resolve the kinetic energy, e^, of the rf out-coupled 




FIG. 2: Extracting the 3D momentum distribution. 

a A time-of-flight absorption image (~ 145/^mx 145/zm) of 
atoms that have been transferred into a third spin state is 
taken after applying an rf pulse to a Fermi gas on the BEC 
side of the Feshbach resonance, b After performing quadrant 
averaging we use an inverse Abel transform to reconstruct 
the 3D momentum distribution. For this particular example, 
a 2D slice at the center reveals a shell-like structure for the 
momentum distribution of the out-coupled atoms. 



atoms we apply an rf pulse that is short compared to the 
trap period. We then immediately turn off the trap, let 
the gas ballistically expand, and measure the velocity dis- 
tribution using state-selective time-of-flight absorption 
imaging, see Fig. [2l Assuming a symmetric momentum 
distribution, we extract the 3D momentum distribution 
of the out-coupled atoms from the 2D image by perform- 
ing an inverse Abel transform. 

We first consider the case of an ideal Fermi gas. To cre- 
ate a very weakly interacting gas we adiabatically ramp 
the magnetic field to the zero crossing of the Feshbach 
resonance. In Fig. [3K, we plot the intensity, which is 
proportional to the number of atoms transferred into 
the third spin state, as a function of the original single- 
particle energy E s and wave vector k. The data are 
obtained by varying the rf frequency and counting the 
out-coupled atoms as a function of their momenta. We 
define zero energy to be the energy of a non-interacting 
atom at rest in the initial spin state. The intensity map 
for a non-interacting Fermi gas is expected to show delta 
function peaks at E s — tu- The white asterisks mark the 
centers of the intensity at each value of k as determined 
from Gaussian fits; these show good agreement with the 
expected dispersion (black line). The rms width in E s 
of the measured spectrum in Fig. [3^i is 2.1 kHz and is 
due to an energy resolution that comes from the rf pulse 
duration. 

To create a strongly interacting Fermi gas we adiabati- 
cally ramp the magnetic field to the peak of the Feshbach 
resonance where the s-wave scattering length a diverges 
and the dimcnsionless interaction parameter l/k F a = 0. 
Here k F is the Fermi wave vector that corresponds to 
the peak density of the original weakly interacting gas. 
Previous measurements have shown that after the ramp 
to l/k° F a = 0, our Fermi gas will be at (0.9 ± 0.1) ■ T c 
for the superfluid state [2j. With photoemission spec- 
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FIG. 3: Single-particle excitation spectra obtained us- 
ing photoemission spectroscopy for ultracold atoms. 

Plotted are intensity maps (independently scaled for each 
plot) of the number of atoms out coupled to a weakly- 
interacting spin state as a function of the single-particle en- 
ergy E s and wave vector k. The black lines are the expected 
dispersion curve for an ideal Fermi gas. The white points 
(*) mark the center of each fixed energy distribution curve, 
a Data for a very weakly-interacting Fermi gas. The Fermi 
wave vector k% is 8.6 ± 0.3 /im" 1 . b Data for a strongly in- 
teracting Fermi gas l/k° F a = and T « T c . The white line is 
a fit of the centers to a BCS-like dispersion, c Data for a gas 
on the BEC side of the resonance where l/k%a « 1 and the 
measured two-body binding energy is h ■ (25 ± 2 kHz). We 
attribute the upper feature to unpaired atoms and the lower 
feature to molecules. The white line is a fit to the centers 
using a quadratic dispersion. 



troscopy on the strongly interacting gas we extract the 
intensity map shown in Fig. [3p. The interactions lower 
the overall energy and flatten the dispersion curve. In 
addition, the energy width is broadened well beyond our 
energy resolution. 

There is now a wide consensus that interpretation of 
previous rf spectroscopy measurements [H, [ljl 17 1 in 
terms of a pairing gap is a difficult problem that is 
still unsolved theoretically [13] ■ The photoemission spec- 
troscopy technique presented here directly measures the 
occupied single-particle density of states and is there- 
fore well-suited for measuring pairing gaps. In BCS the- 
ory the gap vanishes at T c ; however, in the BCS-BEC 
crossover a pseudogap due to preformed pairs is pre- 
dicted to exist above T r 
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Single-particle energy (kHz) 
FIG. 4: Energy distribution curves for a strongly inter- 
acting Fermi gas. We plot the intensity for selected values 
of k. Each curve is the average for seven neighboring values 
of k in Fig. [3p. The data have been smoothed with a 1.5 kHz 
wide filter. 



almost exactly to a "BCS-like" dispersion curve where 
the BCS gap was replaced by the pseudogap 



11(]. As 



111 . [28[ . Perali et al. calculated 



the spectral function for a homogeneous Fermi gas near 
T c and found that the peaks of the spectral function fit 



a first step to analyzing our data, we fit the centers of 
the intensity at each value of k to this BCS-like disper- 
sion curve [11[, E s = // — y/ (e^ — [i') 2 + A 2 . Here, the 
fitting parameters are the renormalized chemical poten- 
tial \J and the pseudogap A. The best fit, shown as the 
white curve in Fig. [3p, gives // = h - (12.6 ±0.7 kHz) and 
A = /i-(9.5±0.6 kHz). In Fig. 14ofRef. [HI , Perali et al 
also plot an example of predicted spectral functions for 
a few values of wave vector k. To facilitate comparison 
with theory, in Fig. [4] we show measured energy distri- 
bution curves (EDCs) for select values of k. It should be 
noted that in all trapped gas experiments, the density is 
inhomogeneous and the pairing gap will depend on the 
local Fermi energy. Therefore, our data should eventu- 
ally be compared with a theory that includes the effect 
of the trapping potential through, for example, a local 
density approximation. Finally, we note that we have 
performed photoemission spectroscopy for a gas cooled 
below T c (initial T/T F = 0.10) and found that the data 
is qualitatively very similar to that in Fig. [3p. 

Far on the BEC side of the resonance, for l/k F a 3> 1, 
the pairing gap eventually becomes a two-body binding 
rather than a many-body effect that depends on the local 
Fermi energy. We measure the excitation spectrum for 
the Fermi gas at l/k F a w 1 where the molecule binding 
energy measured for a low density gas is h ■ (25 ± 2 kHz) . 
We observe two prominent features, see Fig. [3p. The first 
feature is narrow in energy, starts at zero energy, and 
follows the quadratic dispersion expected for free atoms 
(black line). We attribute this feature to unpaired atoms, 
which may be out of chemical equilibrium with the pairs. 
The second feature is very broad in energy, is shifted to 
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lower energy, and trends towards lower energy for increas- 
ing k. This feature we attribute to atoms in the paired 
state. An excitation gap separating the two features is 
evident in the data. We fit the centers of the molecule 
feature to a quadratic dispersion (white line) with the 
free parameters being the energy offset and an effective 
mass m* . In the BEC limit, where one has tightly bound 
molecules, we would expect the energy offset to be the 
molecule binding energy, which equals 2 A, and the effec- 
tive mass to be —m. This negative effective mass reflects 
the fact that out coupling an atom at high k leaves be- 
hind an excitation in the form of an unpaired atom. The 
best fit to the data gives an energy offset of 28 kHz and 
m* = —1.25 to. 

The large energy width seen in Fig. [3J: is likely due to 
center-of-mass motion of the pairs. For comparison with 
the data, we have performed a simple Monte Carlo sim- 
ulation assuming a thermal distribution for the center- 
of-mass motion and using the predicted distribution of 
relative kinetic energy for rf dissociation of weakly bound 
molecules [2Cj . We assume the pairs are in thermal equi- 
librium with the unpaired atoms and use the measured 
temperature of the rf out-coupled atoms corresponding 
to the upper feature in Fig. [3b. Assuming a molecule 
binding energy of h ■ (25 kHz), the calculation gives the 
intensity map shown in Fig. (5^. 

The occupied density of states is obtained by summing 
the data in Fig. [3]over all k, see Fig. [5j>d. For the nearly 
ideal Fermi gas data, Fig. [5b, we find good agreement 
with the expected density of states for a Fermi gas at 
T = 0.18 Tp in a harmonic trap (red curve). For the 
strongly interacting gas, Fig. [5k, the occupied density of 
states becomes wider in energy and the peak shifts to- 
wards lower energies by an amount comparable to Ep. 
From previous measurements of the in-trap size of the 
cloud [29| we estimate the Fermi energy of the strongly 
interacting gas to be h ■ (12.4 ± 0.7 kHz) (dashed line). 
For the BEC side of the resonance, Fig. 5d, a pairing 
gap between bound pairs and free atoms is readily ap- 
parent. The red curve is the expected density of states 
determined from the simulation of a thermal distribu- 
tion of weakly bound molecules (Fig. [5^)- The only free 
parameter in the simulation is an overall scaling factor. 

In this work, we have used photoemission spectroscopy, 
accomplished by momentum resolving the out-coupled 
atoms in rf spectroscopy, to probe the occupied single- 
particle density of states and energy dispersion through 
the BCS-BEC crossover. In the future, it may be possi- 
ble to use spatially resolved photoemission spectroscopy 
to probe the local pairing gap. Another extension of this 
work will be to study the BCS-BEC crossover as a func- 
tion of temperature and/or unbalanced spin population. 
Photoemission spectroscopy for ultracold atoms is a pow- 
erful and conceptually simple probe of strongly correlated 
atom gases that could be applied to many other atom gas 
systems. In the studies presented here, the atoms are 
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FIG. 5: The occupied single-particle density of states. 

a A calculated intensity map for a thermal distribution of 
weakly bound molecules is shown. The white line is the fit 
to the data shown in Fig. [3F b The density of states for 
a weakly-interacting Fermi gas (black line) agrees well with 
a fit (red curve) for a Fermi gas in a harmonic trap. The 
fit, whose only free parameter is the amplitude, includes our 
measurement resolution. The dashed black vertical line shows 
Ef- c The density of states taken at the peak of the Feshbach 
resonance is shifted to much lower energy, d The density of 
states on the BEC side of the resonance has two features; a 
narrow peak due to unpaired atoms and a broader feature due 
to molecules. The red curve is the expected density of states 
from the simulation shown in a. 

interacting via isotropic s-wave interactions and there- 
fore considering different directions of the out-coupled 
atoms' momenta was not necessary. However, like angle- 
resolved photoemission spectroscopy (ARPES) for solids, 
this technique could also be applied to non-isotropic sys- 
tems such as atoms in an optical lattice, low dimensional 
systems, or higher partial wave pairing of atoms [30j ] . 

Methods 

We evaporatively cool an equal mixture of 40 K atoms 
in the I/,™/) = 9/2, —7/2) and \f,m f ) = 9/2, —9/2) 
states, where / and m/ give the hyperfine level in the 
ground-state manifold. We cool the ga s to in an optical 
dipole trap as described previously [291 ]. The frequencies 
for the cylindrically symmetric trap are f r = 233 Hz and 
f z = 19 Hz. The time-of- flight imaging uses a beam 
that propagates along the z direction. At the end of 
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the evaporation, we adiabatically increase the interaction 
strength by lowering the magnetic field, at a rate of 0.1 
G/ms, to a value near the Feshbach resonance located 
at 202.10 ± 0.07 G [2]. The magnetic field values for 
the data in Fig. [3k,-c are 208.43, 202.10, and 201.51 G, 
respectively. 

For photoemission spectroscopy, we apply an rf pulse 
with a Gaussian amplitude envelope with a 1/e 2 width 
of 240 /is, to transfer atoms from the 1 9/2, —7/2} state to 
the |9/2, —5/2} state. The rf frequency is approximately 
47 MHz. Atoms in the 1 9/2, -5/2} state have a two-body 
s-wave scattering length that is 130 Bohr radii with the 
|9/2, -7/2} state and 250 Bohr radii with the |9/2, -9/2} 
state. Immediately after the rf pulse we turn off the opti- 
cal trap and let the atoms expand for 3 to 6.5 ms before 
taking a resonant absorption image of the 1 9/2 , — 5/2} 
atoms. Typically four images are averaged for each rf 
frequency. For the weakly-interacting Fermi gas, the rf 
power was chosen to achieve maximum transfer. For data 
on resonance and on the BEC-side of resonance, no more 
than 30% of the atoms were transferred into the third 
spin state. At the higher rf frequencies, we increased the 
signal by increasing the rf power and then scaled the data 
to give the appropriate intensity. 

In Eqn. [TJ we take 4> to be the rf resonance energy 
measured for a weakly interacting gas. There is an un- 
certainty of ± 1 kHz in <j) and therefore also in the zero 
of E s . The measured wave vector k has an uncertainty of 
approximately 5% from uncertainty in the magnification 
of the imaging system. 
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